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The NuTeV collaboration reported anomalously large weak mixing angle sin^Ow in comparison 
with the standard model prediction. Neutrino and antineutrino charged- and neutral-current events 
are analyzed for extracting sin^Ow- Although the Paschos-Wolfenstein relation is not directly used 
in the analysis, it plays an important role in the determination. Noting that the target nucleus, iron, 
is not an isoscalar nucleus, we derive a leading-order expression for a modified Paschos-Wolfenstein 
relation for nuclei, which may have neutron excess. Then, using charge and baryon-number conser- 
vations for nuclei, we discuss a nuclear correction in the sin^Ow determination. It is noteworthy 
that nuclear modifications are different between valence up- and down-quark distributions. We show 
this difference efi'ect on the NuTeV sin^Ow deviation. 

PACS numbers; 13.15.-|-g, 13.60.Hb, 24.85.-l-p 



I. INTRODUCTION 

The weak mixing angle siri^Ow is one of the impor- 
tant quantities in the standard model. In the on-shcU 
scheme, it is related to the W and Z masses by sin'^dw = 
1 — m^/m^. Collider experiments provide accurate val- 
ues for these masses and the angle. According to a global 
analysis [|, it is sin'^e^-"^"" = 0.2227 ± 0.0004 by ex- 
cluding neutrino-nucleus scattering data. 

The NuTeV collaboration (G. P. Zeller, K. S. McFar- 
land et. al.) reported recently |^ that the mixing angle 
should be significantly larger: 

sin^0w = 0.2277 ± 0.0013 (stat) ± 0.0009 (syst) , (1.1) 

by using their neutrino and antineutrino scattering data. 
For extracting sin?9w, it is known that the Paschos- 
Wolfenstein (PW) relation |] 
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is useful because uncertainties from charm production 
and possible nuclear effects are much reduced. Here, crj^/^ 
and cr'oQ are the deep inelastic cross sections for neutral- 
current (NC) and charged-current (CC) neutrino inter- 
actions with the nucleon. The factor p, which is the rel- 
ative strength between the neutral and charged currents, 
is taken as one. The NuTeV collaboration measured 
charged and neutral current ratios, R^, = cJ^^/tTp^ and 
= '^i^/'^cc ^^'^ then a Monte Carlo simulation is 
used for relating the data to siri^Ow- Fitting these ratios 
simultaneously, they end up using the PW-like relation 
although it is not directly employed in the analysis B . In 
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this sense, it is mentioned in Ref. ||] that "the NuTeV re- 
sult derives sin^dw from the Paschos-Wolfenstein". The 
result suggests that the left-handed neutral current cou- 
pling should be smaller than expected. If it is confirmed, 
it should lead to a new physics finding |^ . The situation 
is summarized in the paper by S. Davidson et. al. 

On the other hand, there are suggestions from a con- 
servative point of view. Miller and Thomas commented 
[Q that the anomalous result could be explained by the 
shadowing difference between neutral and charged cur- 
rent reactions by using a vector meson dominance (VMD) 
model. The NuTeV collaboration replied to their com- 
ments |5] that the explanation is not favored because the 
shadowing effects are subtracted out in the PW relation. 
Furthermore, the model cannot explain observed R^, and 
i?p ratios, which are smaller than expected, and also the 
VMD model does not have proper dependence. How- 
ever, their discussion is refuted by Melnitchouk and 
Thomas in Ref. Nuclear corrections are also dis- 

cussed by Kovalenko, Schmidt, and Yang by noting 
nuclear modifications of -F2- However, one should note 
that such nuclear effects were taken into account in the 
NuTeV analysis in a slightly different way | |To[ |. 

It is not the purpose of this paper to examine the de- 
tails of these previous studies. We rather try to address 
ourselves to the extraction of sin^0w from nuclear data 
in a model independent way as much as possible by re- 
sorting to charge and baryon-number conservations. Be- 
cause the NuTeV target is mainly the iron nucleus, nu- 
clear corrections should be carefully taken into account 
for a precise determination of siri^Ow- In this paper, we 
derive a modified PW relation for general nuclear targets. 
Then, we discuss a possible nuclear modification factor 
which could change the extracted sin^dw value. 

This paper consists of the following. First, nuclear cor- 
rections of the PW relation are discussed in Sec.|l|. Then, 
possible effects on the extraction of sin^Ow are explained 
in Sec. [II. The results are summarized in Sec.^ 
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II. MODIFIED PASCHOS-WOLFENSTEIN 
RELATION FOR NUCLEI 

The PW relation was derived for the isoscalar nucleon; 
however, the used NuTeV target is mainly iron, which is 
not an isoscalar nucleus. The neutron excess may cause 
unexpected nuclear corrections, which should be care- 
fully investigated. In this section, we derive a leading- 
order (LO) PW expression for general nuclei in a model 
independent way. 

First, neutrino- and antineutrino- nucleus charged cur- 
rent cross sections are given in the LO by |]lT| 
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Here, ctq is defined as ctq = Gps/ir, by neglecting the 
factor Q'^ /M^ from the propagator, with the Fermi cou- 
pling constant Gp and the center-of-mass squared energy 
s. The variables x and y are defined by the momentum 
transfer square (= — g^), the energy transfer q*^, and 
the nucleon mass M as a; = Q^/(2Mq°) and y — q^/Ei, 
or jEy. Nuclear quark and antiquark distributions 
are denoted q^{x) and q^{x), respectively. They de- 
pend also on Q^; however, the explicit dependence is 
abbreviated in the parton distribution functions (PDFs) 
throughout this paper. Next, the neutral current cross 
section for neutrino scattering is given by 

"'^NC ^ r r„,2 , „,2 ^1 „.\2\t„,A, 
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+ {4 -f- 4(1 - yfUd^ix) + s^ix)} ] , (2.2) 



where left- and right-handed couplings are expressed by 
the weak mixing angle as ul = 1/2 — (2/3) sin^Bw, 
UB, = -{2/3) sin'^ew, di = -1/2+ {1/3) sin^Ow, and 
dn = (1/3) sin^Bw- The cross section for antineutrino 
scattering is obtained by exchanging the left- and right- 
handed couplings. 

Using these equations, we obtain a nuclear PW relation 

as 
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Here, the valence quark distributions are defined by q. 
qA_qA^ The name , valence quark, is conventionally used 
for valence up- and down-quark distributions, but we use 



the same nomenclature for strange and charm distribu- 
tions by applying the same definitions, = — and 
= — c^. Of course, there is no net strangeness and 
charm in ordinary nuclei, so that "valence" strange and 
charm distributions should satisfy j dxs^{x) ~ and 
J dxc^{x) = 0. However, these equations do not mean 
that the distributions themselves vanish: s^{x) ^ and 
c^{x) + 0. 

The valence-quark distributions and d^ are ex- 
pressed by the weight functions Wu^ and Wd^ at any Q^: 



u^{x) = Wu^{x,A,Z) 



Z Uy{x) + N dy{x) 
A ' 



Mr \ (A r^-, Z dy{x) + N Uy{x) 

dy(x) =Wd^(x,A,Z) J , (2.4) 

where Uy and dy are the distributions in the proton, Z 
is the atomic number, N is the neutron number, and A 
is the mass number of a nucleus. The weight functions 
are defined at fixed (=1 GeV^) in Ref. ||l2[; however, 
they are used at any throughout this paper. The 
explicit dependence is abbreviated as for the PDFs. 
Although these equations are originally motivated by the 
isospin symmetry in nuclei for the virtual Wu„ = = 1 
case, we do not rely on such an assumption. This is 
because nuclear modifications, including isospin violation 
Ip^ , |l^ and nuclear charge-symmetry breaking ||l^ , could 
be taken into account by the weight functions in any case. 
Therefore, the expressions are given without losing any 
generality. 

Next, we define neutron excess constant in and a re- 
lated function Sn^x) by 

N-Z ^ Uy{x)-dy{x) 

en{x) ^ En — (2.5) 
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Then, a difference between the weight functions is defined 

by 

/ X _ Wrf„ {x, A, Z) - (x. A, Z) 
'^^''> wdAx.AZ) + WuA^,A,Z)- ^^■''> 

The function Sy depends also on A, Z, and Q^, but these 
factors are abbrev iate d in writing £„. Substituting Eqs. 
( ^.4| ), (2.5), an d (|2.6|) together with the coupling con- 
stants into Eq. (|2.3|), we obtain 
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Here, Ss and Ec are defined by Eg = s^/[wy {uy + dy)] 
and Ec = Cy/[wy {uy + dy)] with Wy = (wd„ -I- Wu^)/2. 
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We would like to stress that the LO expression in Eq. 
(2.7) has been derived without using any model depen- 
dent factor and and any serious approximation. 

The strange quark (s^) effects are investigated in Ref. 
iQ, and the neutron-excess effects, namely the e„ terms, 
are taken into account in the NuTeV analysis . In 
addition to these corrections, we notice in Eq. (2.7) that 
another correction factor e„ contributes to the deviation 
from the PW relation 1/2 — sin^Ow due to the difference 
between the valence up- and down- quark modifications in 
a nucleus. This fact needs to be clarified. There are 
two restrictions on the valence-quark distributions in a 
nucleus. One is the baryon number conservation [^ , 
and the other is the charge conservation Nuclear 
baryon number and charge have to be A and Z, and 
they are expressed in the parton model as 

A^l dxAY,\{q^ -q^)^l dx^{u^ + d^) , 
Z^j dxAj2eq{q^~q^)^J dx^{2ut-d^), (2.8) 



where A is multiplied in the integrands because the nu- 
clear parton distributions are defined by those per nu- 
cleon, and the relation J dx Sv{x) = Jdxcv{x) = is 
used in obtaini ng t he right-hand sides. Substituting Eq. 
(p^ into Eq. (O), we obtain 



correction of Sy in Eq. (2.7), we obtain 
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where Aw^ is defined by Aw^, — Wy ~ \. Although it is 
not s traig htforward to determine £v{x) from Eqs. ( |2.9| ) 
and (2.10), it indicates that ey{x) is finite. 

In this way, we find that nuclear modifications are, 
in general, different between the valence up- and down- 
quark distributions because of the baryon number and 
charge conservations. It gives rise to the factor £„. This 
kind of detailed nuclear effects cannot be accounted sim- 
ply by investigating electron and muon scattering data 
and by fitting charged current cross-section data for the 
same target . Because the physics associated with the 
Ey factor is missing in the NuTeV analysis, it may cause 
a significant effect on the sin'^dw determination. 



III. EFFECTS ON sin^6lw DETERMINATION 



The angle sin^Ow can be extracted by using Eq. (2.7) 



together with the experimental data of i?^ . In order to 
find whether or not the e„ corrections could explain the 
deviation from the standard model value for sin'^Bw, we 
approximate the relation by considering that the correc- 
tions are small ( Sy <^ 1). Retaining only the leading 
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The higher-order and other corrections O(e^), 0(en), 
0(es), and 0{ec) are not explicitly written. As far as 
present neutrino data suggest, the "valence" strange dis- 
tribution should be small, and the measurements indicate 
that such a correction increases the NuTeV sin^Ow de- 
viation according to Ref. ff4|. Therefore, at least at this 



stage, the finite s^ and cj^ distributions effects, 0{£s) 
and 0{£c), are not the favorable explanation. Although 
accurate measurements may clarify the details of the dis- 
tributions s^ and Cy in future, they are not discussed in 
the following. The neutron-excess effects 0(£„) are in- 
cluded in the NuTeV analysis |jl^, so that they are not 
the source of the sin^Ow deviation. 

As mentioned in Sec. | and Ref. |5|, the NuTeV 
sin^Ow is "derived" from the PW-like relation indirectly. 
It is obvious from Eq. (3.1) that there is an £„-type 
correction to sin^Ow, and it may explain, at least par- 
tially, the deviation from the standard model. However, 
the correction is essentially unknown at this stage in the 
sense that there is no significant data to find the differ- 
ence between valence up- and down-quark modifications 
{sy). In order to investigate whether or not the £y cor- 
rection is large enough, we should inevitably use some 
prescription for describing the Sy factor. In the follow- 
ing, we introduce two different descriptions as examples. 



1. A prescription for the conservations 

It is not str aightf orwar d to find a solution £y{x) to sat- 
isfy Eqs. and ( pTO| ). For an approximate estimate, 
the higher-order corrections £y£n a re n eglected in these 
equations. Then, substituting Eq. (2.9) into Eq. (2.10), 
we obtain 



ey{x) 



Uy(x) — dy{x) /S.Wy{x) 
Uy{x) + dy{x) Wy{x) 



(3.2) 



by considering a special case that the integrand vanishes. 
Of course, this is not a unique solution, but this estimate 
should be able to provide information about the magni- 
tude of the correction. 

2. A analysis of nuclear PDFs 

Global analysis results could be used for calculating 
£y{x). A x^ analysis for determining nuclear parton dis- 
tribution functions (NPDFs) is reported in Ref. [|l2j, and 
obtained distributions can be calculated by using subrou- 
tines at the web site in Ref. |l^ . Before using the NPDF 
code, we would like to remind the reader what has been 
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done for the valence distributions. At Q^~l GeV'^, the 
weight functions are assumed to be 
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where i denotes Uy or dy, and 

^u^j by J and d>y are the 
parameters to be determined by the fit. Because it 
is the first analysis attempt for nuclei, the parameter 
number has to be reduced as many as possible. It is the 
reason why the common parameters are chosen for /3^, , h^, 
Cy, and dy. However, in order to satisfy both the charge 
and baryon-number conservations, at least one parameter 
should be different. Because the parameters and a^^ 
are determined so as to satisfy these conservations, this 
description of e«(x) surely satisfies Eqs. (^) and (2.10). 
However, one should note that this e^(S) may not be 
valid if it has more complicated x dependence than the 
one calculated from Eq. (3.3). 



A. Numerical results 

In the first description, is evaluated numerically by 
using Eq. (|1|) with the NPDF code in Ref. ^ for 
calculating Uy, dy, Wu^, and Wd^ at given Q^. In the sec- 
ond one, Ey is calculated by the definition in Eq. ( |2.6| ) 
with the weight functions Wu^ and ui^^, which are nu- 
merically calculated by the NPDF code. In the NuTeV 
measurements |] , averages of the kinematical variables 
are given by < £; >=120 and 112 GeV, < >=25.6 
and 15.4 GeV^, < a; >=0.22 and 0.18 for neutrinos and 
antineutrinos, respectively. The ranges of x and are 
0.01 <x < 0.75 and 1 < < 140 QeV^. Although the 
v and i> incident energies are different, we use the aver- 
aged value < E >=116 GeV for connecting x, y, and Q^: 
y = Q^/{2Mx < E >). 

In Fig. |l|, the value is fixed at = 1 GeV^, and 
the correction term in Eq. (^^) is evaluated as a function 
of X. The neutron excess e„ = 4/56 and siri^Ow = 0.2227 
are used in the calculations. The solid {e^v^) and dashed 
[ey ) curves are obtained by the first and second de- 
scriptions, respectively. The dotted line indicates the 
NuTeV sin'^Ow deviation (0.2277-0.2227=0.0050) just 
for a comparison. Both curves increase rapidly as x be- 
comes larger, and this is mostly a kinematical effect due 
to the factor 1/[1 — (1 — yY]. Another effect at large x 
comes from the Fermi-motion-like factor 1/(1 — 2;)'^" in 
Eq. (3.3). The £y correction term is very sm all w ith 
the following reason. The x dependence in Eq. (3.3) in- 

dicates that Sy is almost independent of x except for 
large x. The first terms in Eqs. (2.£) and (2.1C) are 
valence-quark distributions multiplied by Awt, , and these 
integrals almost vanish. The obtained Sy is roughly pro- 
portional to these small integral values. One should note 
that this result could be artificially small due to the weak 
X dependence assumption. On the other hand, the mag- 
nitude of el^^ is large in general and it has completely 



different x dependence. Because e]P is directly propor- 
tional to the valence-quark modification (Awi,) in Eq. 
( |3.2D , the function changes sign at x « 0.3, which is the 
transition point from antishadowing to the EMC-effect 
region. In the revised analysis ||l^, |l9j , the valence modi- 
fications are slightly different; however, the essential fea- 
tures are the same. According to the first description, the 
function becomes comparable magnitude to the NuTeV 
deviation. 
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FIG. 1: The £v correction term is eval uate d at Q^=l GeV^. 

The solid curve is calculated by Eq.^ ( p.2[) , and the dashed 

curve is obtained by the fit code The NuTeV de- 

viation 0.005 is shown by the dotted line just for comparison. 
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FIG. 2: The correction term is evaluated at Q'^=20 GeV^. 



Because the average is much larger than 1 GeV^ 
in the NuTeV experiment, the corrections are also cal- 
culated at (5^=20 GeV^ by noting the kinematical limit 
y < \. The results are shown in Fig. ^. In compari- 
son with the Q^=l GeV^ results, the effects are much 
suppressed. This is again due to the kinematical factor 
1/[1 — (1 — y)^]. For example, this factor is 55 for x=0.5 
and Q^=l GeV^; however, it becomes 3.0 for a;=0.5 and 
Q^—2Q GeV^. This is the reason why both distributions 
become smaller. Although the el, effect is too small to 
explain the deviation at Q^=20 GeV^, the ei^"* is still 
comparable magnitude. 

In oder to investigate these effects on the NuTeV 
sin?6w, the analysis should be done with the Monte 
Carlo code with the experimental data. However, in or- 
der to show the order of magnitude, we simply average 
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the obtained curves over the x range (Ax) from 0.01 or 
Xmin ~ Q'^/{2M < E >) to 0.75. Because the data are 
centered at about x — 0.2, this kind of simple average 
could overestimate the effects coming from the large x 
region. The calculated results are shown by the solid 
curves in Fig. ^. As aheady found in Figs. |l| and ^, the 
effects are very large at small and become smaller as 
becomes larger. 




1 10 100 



Q' (GeV') 

FIG. 3: The solid curves indicate the corrections averaged 
over the x range. The dashed ones are obtained by taking the 
NuTeV kinematics into account. 

If the simple x average is taken, the effects look large. 
On the other hand, there is a method to take the NuTeV 
kinematics into account |p2| by using the functionals in 
Fig. 1 of Ref. jllj. Although the physics motivation 
is completely different, the present Sy distribution could 
be simulated by the NuTeV distributions uf, — c?" and 
— u". It is interesting that their "isospin- violating 
distributions" could effectively contain the present nu- 
clear effect. If such a correspondence is made, the effect 



on siri^dw is calculated by using the NuTeV function- 
als 1^, 1^. The results are shown by the dashed curves 
in Fig.nH We find that the effects become significantly 
smaller due to the lack of large x events in the NuTeV 
experiment. Therefore, as far as the considered two de- 
scriptions are concerned, the nuclear modifications are 
not large enough to explain the whole NuTeV deviation. 

It is, however, still too early to conclude that the 
present mechanism is completely excluded, because the 
nuclear modification difference between and is not 
measured at all. It may not be possible to find this nu- 
clear effect until a neutrino factory ||l^, ^ or the NuMI 
project is realized. On the other hand, if these 
facilities are built, the "nucleon" cross sections (hence 
sin^dw) could be measured with proton and deuteron 
targets with minor nuclear corrections. The NuTeV also 
reported that Ri, and i?p are unexpectedly smaller 
We are also investigating this issue, and hopefully it will 
be clarified in the near future. 



IV. SUMMARY 



We have derived a modified Paschos-Wolfenstein rela- 
tion for nuclei. Using this relation, we investigated the 
possibility that the NuTeV sin^Ow deviation could be 
explained by the nuclear parton distributions in iron. In 
particular, we pointed out that nuclear modifications are 
different between valence up- and down-quark distribu- 
tions. The difference partially explains the NuTeV de- 
viation although it may not be large enough to explain 
the whole deviation. Because such modifications are not 
measured, it is important to investigate them in future. 
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